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Abstract 

' In the present paper which follows our previous paper "Mathematical 

models for passive imaging I: general background" , we discuss the case of 
^sO ! surface waves in a medium which is stratified near its boundary at some 

scale comparable to the wave length. We discuss how the propagation of 
such waves is governed by effective Hamiltonians on the boundary. The 
results are certainly not new, but we have been unable to find a precise 
i -pH | reference. They are very close to results in adiabatic theory. 

Introduction 

>< : This paper is strongly related to the first part . We will be more specific and 

discuss the case of surface waves which are used in seismology in order to image 
the earth crust. 

We consider a medium X with boundary dX and assume that the medium is 
stratified near dX. We will discuss how the linear propagation of waves located 
near the dX is determined by an effective Hamiltonian on dX. It is interesting 
enough to remark that this Hamiltonian is no more a differential operator, but 
only a "pseudo-differential operator" (a \l/DO) with a non-trivial dispersion 
relation (principal symbol). This situation is well known in physics as giving 
birth to some kind of wave guides. 

A typical motivating situation is that of seismic surface waves propagating 
along the earth crust: it is well known that, at "macroscopic scales", the earth 
crust is horizontally stratified such giving birth to the so-called surface seismic 
waves which admit some more technical names like "Rayleigh" or "Love" waves 
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(see P). They play a crucial role in the passive imaging of the earth crust. 
Moreover, they are the most dangerous in case of an earthquake. 

From the general methods described in our paper the field-field correlation 
allows to recover the dispersion relation of these surface waves from which we 
want to recover the transverse (vertical) structure of the crust. After describing 
a general result which is close in spirit to the adiabatic Theorem in Quantum 
Mechanics, we will discuss briefly what kind of inverse spectral problem we need 
to solve at the end and we solve it under a physically realistic monotonicity 
assumption. 



1 The general setting 

We will work locally in X = {(x, z) G ]R d_1 x R \ z < 0}. We will consider the 
very simple case of an acoustic wave equation near the origine of X: 

j u tt - div(n gradw) = , , 

\u(x,0)=0 [) 

with 

n(x, z) = iV(x, z, -) 

and iV(x, z, Z) : x R_ x K._ — > IR + a non negative function which is inde- 
pendent of Z for Z < Zq < 0. There are 2 possible assumptions on the regularity 
of N: 

• Assumption 1 N is smooth 

• Assumption 2 N is smooth as a function o/(x, z) with values into L°°(M_) 
and N(x, z, Z) = JV(x, 0, Z) + 0(z°°) 

In other words, the medium X admits a quite irregular behaviour in the 
vertical direction in a very small horizontal band B £ C X of width eZq. We plan 
to see that Equation (JTJ) admits, as e — > 0, asymptotic solutions of frequency of 
order e^ 1 located in B e . Moreover these solutions are determined by solving an 
effective pseudo-differential equation on the boundary dX = x {0}. 

Assumption 3 We will assume that 

< inf N(x, O, Z) < N(x, 0, -oo) . 

Z<0 

Physically, it means that the propagation speed at some points very close to the 
boundary is smaller than inside the medium. This kind of assumption is usually 
satisfied in seismology where the speed of elastic waves into surface sediments 
layers is smaller than the speed inside the rocks below the sediments. 
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2 The main result 



2.1 A Sturm-Liouville operator 

Let us consider, for each (x, £) e T*dX, the self-adjoint differential operator L Xj g 
on the half line Z < 0, with Dirichlet boundary condition at Z = 0, defined by: 

^ := -fi{N(x, 0, Z) — ) + iV(x, 0, Z)|£| 2 ?, (2) 

In case of Assumption |2J £ x .£ is defined in terms of quadratic forms. 

The spectrum of L x ,t consists of a finite discrete spectrum and a continuous 
spectrum [iV(x, 0, Z )|£| 2 , +oo[. Under Assumption 03 L x ,c admits, for £ large 
enough, a non empty discrete spectrum of simple eigenvalues 

inf A(x, 0, Z) \i\ 2 < Ax(x, f) < • • ■ < A, (x, £)<-■-< A fc (x, £) < A(x, 0, Z ) |£| 2 , 

which depend smoothly of (x, £). In order to see that, we can interpret |£| _2 L Xj £ 
as a semi-classical Schrodinger type operator with an effective Planck constant 
l^l" 1 and a principal symbol 

Px (Z,C) = A(x,0,Z)(C 2 +l) 

which admits a well near (Z, £) = (0, 0). We should however take care of the fact 
that the number k depends on (x, £) and goes to oo as £ does. 

It leads to an interesting bifurcation problem which we will discuss in Sec- 
tion m 



2.2 WKB solutions of the stationnary wave equation 

The goal of this section is to build the effective surface Hamiltonians which de- 
scribe the surface waves. Let us start with the: 

Lemma 1 Let us consider the operator H defined by: 

Hu := — 5 2 div(n gradw) (3) 

acting on functions on X vanishing at z = (Dirichlet boundary conditions). 

Let us choose A(x, £) an eigenvalue of L^ depending smoothly o/(x, £) e U, 
where U is a bounded open set ofT*dX, and y?(x, £, .) a normalized associated 
eigenf unction. There exists: 

• An asymptotic expansion 

oo 

<£ £ = 5> m (x,£,-)e m 

m=0 

with ipo = tp and the y m (x, £, .)'s are smoothly dependent of (x, £) 
values in the domain of L x ^. The tp m 's (m > 1) are unique if they are 
assumed to be orthogonal to ip. 
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• A symbol 

oo 

a B (x,0 = ^a m (x,0^ m 

m=0 

such that we have the following identity of formal power series in e: 

h ($ £ (x, e, ^y<*K>/ £ ) = a £ ( X , o$ £ (x, e, . (4) 



Proof. 



Expanding N by Taylor formula, we get: iV(x, eZ, Z) = Xl^o -^( x > 
with 

1 d l N 

Under Assumption Q for / > 1, the JVj's are smooth and compactly 
supported in Z. Under Assumption |2J the iVj's vanish for I > 1. 
We see that we are reduced to compute 

: = fl-(e*W€>/V(x,^)) , 

with n(x, z) = iV(x, z/e). We find 

m = e *(*IO/e (L Q <p + eL x ip + e 2 L 2 ^) 

with L = L x ,§, 

= -2m(£|V x y?) - i(£|V x n)y9 

L 2 <£> = -(V x n|V x y?) > -nA x tp . 

We start with tp = tp an eigenfunction of L x g with eigenvalue A. 
By induction, we pick all terms in e m and get from Equation (jlj): 

L x £<p m = Ay?m + Om^O + R(<fO, " ' ' , ^m-l) , (5) 

with 

TO— 1 

-R(^p, • • • , Vm-i) = flm-iyi - y Iy<# 

i=l l<m, k+l+j=m 

where is given as Lj with A replaced by N^. 
We first choose the unique a m so that 

a m ^ + R(fo, ••■ , (Pm-l) 

is orthogonal to tpo, then we can solve Equation (J5J). In order to 
show that <£> m lies in the domain of £ x ,£, it is enough to show that 
R(<p , ■ ■ ■ ,<p m -i) is in L 2 (Z). 



□ 
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2.3 A short review on micro-functions 

In order to help the reader, we give here a very short review on micro-functions 
and $DO's. Good references for more details could be [2j or 
Let us first give the main definitions: 

Definition 1 A family of functions (distributions) f £ : X — > C is said to be 
admissible if for any function \ £ C^°(X), there exists some real number s so 
that the Sobolev norms ||x/ e || s are at most of polynomial growth w.r. to e~ l . We 
will denote A(X) the vector space of such families. 

Definition 2 The frequency set or mi crosupport, denoted WF(f e ) , of an admis- 
sible family f £ is the closed subset ofT*X which is defined as follows in canonical 
local coordinates (x,£): 

K£o) i WF(f £ ) if and only if 3 X G C™(X), x {x ) ± 0, 
TAxfM) = O(e^) for £ close to £o . 

Here T £ u is the e— Fourier transform: 

F e u(0 = (2irs)- d/2 J e- i{xl0/e u(x)\dx\ . 

Definition 3 If X is a smooth manifold and U is an open set in T*X, M.(U) 
the space of microf unctions in U, is the quotient 

M(U) := A(X)/{f £ \WF(f £ ) n U = 0} . 

We will denote f £ = 0(e°°) in U the property 

WF(f £ ) n u = . 

Definition 4 If a(x,£,e) is a suitable function, called a symbol, the pseudo- 
differential operator (^DO) A = Op e (a) is defined by: 

Au(x) := {2ne)- d J e i{x - y ^ /£ a(x,£,e)u(y)\dyd£\ . 

The main result is that \l/DO's act on microfunctions: 

Theorem 1 If a(x,£,e) belongs to a suitable class of symbols (for example, all 
derivatives of a are uniformly bounded independently of e ), then 

WF(Op E (a)(f £ ))cWF(f e ) . 

Moreover the action of Op e (a) on Ai(U) depends only on the values of a in 
some neighbourhood of U . 
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2.4 Reformulation of Lemma [T] in terms of micro- functions 

From Lemma we get the 

Theorem 2 Let us fix (xo, £o) £ T*dX and assume that we have a smooth eigen- 
value A(x, £) of defined in some bounded neighbourhood U' o/(x ,£o)- 
its <jroe ?7, an open neighbourhood of (x ,£o)> 50 U C U' and x £ C™{U') 
which is = 1 on U. 
The map: 




defines, by passing to the quotient, a linear map J £ : ,M(£/) — ► A^(V), where 
V:={(x,z;t,()eT*X \ (x;£)eU} , 

which satisfies 

• J £ is microlocally unitary 

• J £ is independent of the choice of % 

• WF (T £ (u)) n V C {z = 0} 

• if J £ = J £ H c f[, where H e s is a ^>DO of full left symbol a e (x, £) and of 
principal symbol A(x, £) 

Proof. - 

The map J e is well defined because <E> e (x, £, Z) is well defined for 
(x, £) G £/'. Let us explain first how J £ is defined. We first can define 
A by 

A £ = Op £ ( X (x,0$ £ (x,£,.)) , 

as a \l/DO whose symbol takes values in L 2 (M~). As such, by The- 
orem HI it is well defined from M(U) to M(U; L 2 (R~)). We now use 
the unitary map E £ : M(U; L 2 (M~)) — »■ A4(V) which is defined by 
E £ : f(x, .) — > e - 2/(x, z/e) and define J £ = D £ o A £ . 

The rapid decay of $ e (x, £, Z) w.r. to Z implies the property 
f£(u)) n V C {z = 0}. 

The unitarity comes from the symbolic calculus: A* o D* o D £ o A £ 
is a scalar ^DO of principal symbol |x( x ? £) | 2 £ _1 |<K X ) -)\ 2 dz 
which is = 1 on U. 

The interlacing property is easily checked on exponentials and 
hence on microfunctions as a consequence of Lemma ^ 

□ 
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2.5 Speeds of propagation 

To any dispersion relation K(x.,£) = u 2 is associated a speed of propagation (the 
so-called "group-speed") defined by 

V := f = 8 ( K/2, . 

For the acoustic equation, we have n(x, z)|£| 2 = uj 2 and hence v = n^/u and 
i) := |v| = y/n. We see the known result that the speed is independent of the 
frequency 

Let us compute the speed for A(x, £) = u 2 : we have 

«9 5 A(x, = 2 (7° iV(x, 0, Z) V (Z) 2 dZ^ i 



and hence 



N(x,0,Z)ip(ZydZ 



and using 



r° 

/ iV(x,0,Z)|£|V(Z) 2 dZ < A = uo< 

J — oo 



we get 

The speed of the surface waves is less than the speed of the body waves. 



v < — < v/ra(x, 0, Z ) . 



3 An inverse spectral problem 

Recovering the vertical structure of earth crust from the dispersion relation of 
surface waves needs to solve the following inverse problem: 

Let us assume that we know the discrete spectrum of L x ^ at fixed 
x for all £'s, can we recover the vertical profile iV(x, 0,Z)? 

Let us show the 

Theorem 3 If iV(x, 0, Z) is a decreasing function of Z from iV(x, 0, Z ) to 
iV(x, 0,0), the asymptotics of the discrete spectra Aj(x, £), 1 < j < fc(x, £) as 
£ — > oo determine the functions iV(x, 0, Z). 

Proof. - 

This is a consequence of Weyl's asymptotics: for E < iV(x, 0, Z ), 
let us introduce 

iV(x,£,E) = #{A J (x,0<^|£| 2 }. 
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We have the following (Weyl type) asymptotics as £ — > oo: 

N(x,£,E)~ (|JY Area({iV(x,0,Z)(l + C 2 ) < E}) . 
It implies that, from the asymptotics of Aj(x, £), we can recover 



o 



E 



V(E) = / i/TTT^r - 1 ^ 



with N(z(E)) = E. We can rewrite 



JV(Z) 



•/Wo 

and introducing 



u 



z'{u) du 



a 



r E , 

Kf(E) := / VW^Ti f(u) du 
Jn 



we get 

-^- 3 ((KoK)f)(E) = 2f(E) 

which allows to recover z'(E)/\fE as the third order derivative of 
K(V)(E). 

□ 

It is a quite interesting problem to try to extend the previous Theorem without 
the assumption that N is decreasing. In this case Weyl formula is no more enough, 
but more refined trace formulae like Gutzwiller trace formula may be usefull to 
give a general answer. 



4 Bifurcations 

What happens when the eigenvalue Afc(x, £) goes to the bottom of the continuous 
spectrum. It is an interesting place when there should be a mode conversion from 
body waves into surface waves or conversely. It could be interesting from the point 
of view of the study of earthquakes whose source is deep: how do they give birth 
to surface waves? It could also be interesting in order to study resonances created 
by surface waves into a bassin. We plan to study these (difficult?) questions in a 
future work. 
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